function draw_feature( fval, classifiers, cid, fid, ths)

hold on;
fval = fval{cid};
classifier = classifiers(cid);
cmap = [0,0,0.6; 1,0,0; 0,0.5,0; 0,0.5,0.6; 1,0,0.6; 1,0.5,0];

for i=1:length(fval.features)
    mu = fval.features{i}.mu(fid);
    sigma = fval.features{i}.sigma(fid);
    x = mu-3*sigma + [0:0.06:6] * sigma;
    plot(x, normpdf(x, mu, sigma), 'Color', cmap(i, :), 'DisplayName',strrep(implode(classifier.classes{i}, ','),'_','\_'));
end

yl = ylim;
for i=1:length(fval.features)
    ypos = yl(1)+(rand(size(fval.features{i}.val, 1), 1)*(-0.13)+0.15*1)*(yl(2)-yl(1));
    plot(fval.features{i}.val(:, fid), ypos, 'o', 'Color', cmap(i, :), 'DisplayName',strrep(implode(classifier.classes{i}, ','),'_','\_'));
end

yl = ylim;
for i=1:length(fval.features)-1
    ypos = yl(1)+[0.1 0.9]*(yl(2)-yl(1));
    plot([ths(i) ths(i)], ypos, 'black');
end

legend('show');

end
